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OJV AD J U8TMENT F0EMULA8. 



EY E. L. DE FOEEST, M. A. 

(GonUnuedfrom page 86.) 

Theee separate tests of good adjustment have been proposed by the pres- 
ent writer, for which see Smithsonian Beport of 1873, p. 333, and a more 
recent essay in pamphlet form j [Interpolation and Adjustment of Series, 
New Haven, 1876). The first of these tests is, that if we denote by v the 
residual error, found by subtracting an adjusted term from the correspoding 
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given one, and denote by e^ the mean error of the given term, then form- 
ing the values of (f-j-ej)^ for all the terms of the series, we ought to have, 
denoting the whole number of terms by N, 

^.{±y = i^.m.^^ (9) 

that is to say, the arithmetical mean of all the values of {v-^-eiY ought to be 
approximately equal to unity, and ought not to diifer from it by more than 
its probable error, which is 



.6745. 



'^|w 



The several values of £j are supposed to be determined from the nature 
of the observations, each given term being the mean result of a number of 
observations. The values of e^ will often vary for different terms of the 
series, not merely in inverse ratio to the square root of the number of ob- 
servations made upon each term, but also as some function of the term itself, 
80 that the mean error of one term cannot always be inferred from that of 
another term solely by comparing the numbers of observations made upon 
each. An illustration of this is furnished by formula (96), in the Smithson- 
ian Beport of 1873, p. 334. 

There is an even chance that the true series, if we had it, would satisfy 
the test of good adjustment just described. If, therefore, we find that our 
adjusted series does satisfy it, we may consider that a good approximation 
to the truth has been reached, and perhaps as good as the nature of the case 
will permit. If the mean of (v-i-SiY falls below the lower limit, that is falls 
short of uniiy by more than its probable error, we must infer that the resid- 
uals V are too small, and that the series probably has not been smoothed out 
quite enough, while if it goes beyond the upper limit, so as to exceed unity 
by more than its probable error, the inference will be that it has been 
smoothed too much. It thus becomes interesting to see whether we cannot 
find some rule to guide us in the choice of a formula which will probably 
make an adjustment satisfying this test. Let us inquire then, what the most 
probable value of the arithmetical mean of (f-r-Si)^ will be, for any given 
adjustment formula used. 

Assuming that all the terms of the given series are of equal weight, let n 
be a number such that the probable error, or deviation from zero, of the 
arithmetical mean of any n consecutive values of v shall be equal to the 
probable error of an adjusted term, so that we have 

s'=-^ (10) 
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The probable error of any one of these n terms will be approximately 

If we denote by v^ the true error, or remainder after subtracting the true 
value of a term from its given or observed value, we shall also have 

. = .6m^(m)). 

Placing these two values of e equal to each other, and giving to n its 
value from (10), we get 

But the most probable ratio of the two sums in the first member will evi- 
dently be unchanged if we extend the summation throughout the whole 
series of iV^ terms, and dividing both numerator and denominator by e|,and 
also by N, we shall have 

Now in the first member the most probable value of the denominator is 
unity, because every value oiv^-^s^ is a true error in a system whose mean 
error is unity, and the square of the mean error is the mean of the squares 
of all the true erorrs. Hence we have, as the most probable value of the 
arithmetical mean of {v-^s.^^ 

^-(~r=;-(4)' (") 

It thus appears that this mean value will most probably be less than unity 
by the square of the error-ratio due to the adjustment formula employed, 
though it may and probably will exceed or fall short of this value in par- 
ticular cases.* We have reached this result by assuming that the weights 
of the given terms are all equal, yet it will hold true approximately if they 
are unequal, for the probable value of v will vary in nearly the same ratio 
as £i does. 

Since the proposed test of good adjustment requires that the mean of 
(t)-i-£j)^ shall not fall short of unity by more than its probable error, which is 



.6745 



*Its probable error, or probable deviation from the most probable value, will be 
Compare Interpolation and Adjustment of Series, p. 10, 
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if the adjusted series is the true one, it follows that we ought to choose an 
adjustment formula whose error-ratio e'-f-e is at any rate small enough to 
satisfy the equation 



(4y=l-..7«V^ 



N 



that is, we should have 
which may be written 



r<iJ(-«'«V|) (12) 

il < -9767 

e = i/N' 
If we wish to adjust a series of 50 terms, we cannot consider that the chan- 
ces are in favor of the test being satisfied, unless we use an adjustment for- 
mula whose error-ratio does not exceed 

.9767 „„, 

so that we should choose a formula of at least as many as 17 terms in Table 
A., or 23 in Table B., or 19 in Table C. On the other hand, formula (12) 
gives 

N^ 2{.e745f-i-(-^y (13) 

that is to say N '—, — r^:. 

Assigning to e'-i-e the successive values .300, .345, .311, 
we find for i^ the values 112, 64, 97, 

showing that while a series of as many as 112 terms can probably be ad- 
justed by the 25-term formula of Table A. so as to meet the test, the long- 
est formula of Table B. would probably succeed only for a series of 64 terms 
or less, and that of Table C. for one of 97 terms or less. 

The conclusion we come to from this investigation is, that it will be 
advisable in most cases to use the longest adjustment formulas given in the 
tables, provided that this can be done without violating the condition that 
any 2m-|-l included terms must not deviate greatly from the form of a se- 
ries of the third or any lower order. 

As the first m and last m terms of the series cannot be reached directly by 
the formula, the series should be graphically extended by m terms at both 
ends, first plotting the observations on paper as ordinates, and then extend- 
ing the curve along what seems to be its probable course, and measuring the 
ordinates of the extended portions. It is not necessary that this extension 
should coincide with what would be the true course of the curve in those parts, 
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The important point is, that the m terms thus added, taken together with 
the m+1 adjacent given terms, should follow a curve whose form is approx- 
imately algebraic and of a degree not higher than the third. As the adjust- 
ing process requires repeated multiplications by the coefficients I, it will be 
well to prepare in advance a table showing the product of each of these co- 
efficients by each of the nine digits. It will also answer every purpose, and 
save labor, if we adjust only alternate terms in the regular way, and then 
fill in the intermediate terms from these by the simple formula for "inter- 
polation into the middle", which may be written 

«o=M9(«i+«-i)— ("s+w-s)] (14) 

Two additional tests of good adjustment have been discussed by the writer, 
one of which depends upon the fortuitous grouping of the + and — signs 
in the series of residual errors v. {Interpolation and Adjudmeni of Series, 
p. 31.) This is much the simplest test of any, its use requiring little labor 
or knowledge of mathematics. It was shown that in the case of a periodic 
series, that is, a series whose first and last terms are consecutive, the moat 
probable number of isolated groups of n like signs, occurring from accidental 
causes, will be 

where the expression which follows the doubtful sign is the probable error. 
This will apply sufficiently well to the case of ordinary or non-periodic se- 
ries, when the method of adjustment here discussed has been used, provided 
that we treat the first and last signs of v as consecutive, so that if they are 
alike, they belong to the same group. Taking w = 1, we find that the most 
probable number of isolated single signs, that is, signs unlike both the ad- 
jacent ones, is 

lN±--^T^V[m), (16) 

and for n = 2, the most probable number of groups of two like signs is 
found to be 

so that the most probable number of signs falling within groups of two is 

lN±'&^x/{'iN) (17) 

Adding (16) and (17) together, we have for the most probable number of 
signs falling within isolated groups of either one or two like signs 

IN+'SJ^x/i^^N), (18) 

which may be written 

IN± .533i/iV, 
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and this is also the expression for the most probable number of signs falling 
in groups of more than two. Hence we have this practical rule, that if a 
series has been well adjusted, the whole number of signs of the residual v 
which fall within groups of only one or two like signs each, will probably 
be about equal to the whole number which fall within groups of more than 
two, and the probable error of either number is 

0.533/iV. 
If the number of signs in groups of one or two should exceed ^N by more 
than this probable error, the inference will be that the series probably has 
not been smoothed out enough, whereas if they fall short of ^Nhy more than 
the probable error, we must presume that the series has been smoothed too 
much. 

It was also shown in the treatise referred to, p. 33, that if in formula (15) 
the most probable number N-r- 2"+^ is found to be less than J for any as- 
sumed value of n, it will show that the odds are against the occurrence of 
any group of so many as n like signs. If we take 

1 ^ N 
4 2"+! 
it gives 

•' = ^ + ^ (") 

that is to say 

«, = 1 + 3.32 log iV. 

Hence we have the rule, that there will be a preponderance of chances against 
the occurrence, from accidental causes, of any group of a number of like signs 
greater than 

1 + 3.32 log JV, 
and if a larger group does occur, it will indicate a probability that the se- 
ries has been smoothed out too much, or that its true law, in that vicinity, 
cannot be fairly represented by an algebraic curve of the third degree, for 
so many consecutive terms as are included by the adjustment formula. 

The writer will take this opportunity to make a remark about the method 
of constructing equations of curves representing annual variations of tem- 
perature, from the monthly means taken as data, referred to at p. 314 of the 
Smithsonian Report of 1871. It was a means of suggesting to my mind a 
general method of interpolation, and has been referred to as the discovery of 
Professor Everett, who published it as such in the Edinburgh New Philo- 
sophical Journal for July 1861, and again in the American Journal of Sci- 
ence and Arts for January 1863. I have since learned that the method had 
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been published some eleven years earlier, by M. Bravais, a French meteor- 
ologist. It may be found at p. 324 of Vol. II. of the Meteorologie which 
forms part of the series of the Voyages de la Commission Scientifique du Nord, 
edited for the French Government by M. Gaimard and others. The vol- 
ume mentioned has no date, but one of the meteorological charts in the ac- 
companying Atlas de Physique bears the date 1850.* The analytical process 
in question, therefore, ought to be designated after its first discoverer, as the 
method of Bravais, rather than as Everett's method. Of course this earlier 
origin of the method increases the probability that Schiaparelli was not un- 
acquainted with it, as suggested by me in Interpolation and Adjustment of 
Series, p. 39. The same property has recently been published by a writer 
in the Astronomisohe Naohriohlen for Nov. 24, 1873, and this is catalogued, 
apparently as a new discovery, in the Jahrbuoh uber die Fortsohritte der 
Mathematik for that year, p. 123. 

An inaccurate remark was made by me at p. 310 of the 8m. Report of 
1871, in describing M. Tchebychef's mode of arranging data as intended 
for making "ordinary interpolations, not from groups, but from single terms 
or ordinates." From a brief allusion to the method in Bertrand's Differ- 
ential CMeulus I had supposed that it was something like what is known as 
the method of normal places. That this was an error, has been shown by 
a recent examination of Tchebychef's original memoir, which is entitled 
8ur V Interpolation dans le ca^s d'un grand nombre de donnees fournies par 
les Observations. The method regards the series of observations as geomet- 
rically represented by ordinates, and the area of the polygon formed by join- 
ing their extremities is supposed to be divided by certain limiting ordinates 
into a number of areas, which are regarded as values of the area Jydx for 
the required curve, between those limits.f These areas are taken as data for 
determining the valuos of the constants in the algebraic equation of the curve, 

2/ = ^0 + ^ix + A^x" + + A«". 

The constant A^is determined from one set of areas, J.^ from another, and 
Ai^ from still another, and so on, the limits of these areas being chosen with 
a view to securing the best values for the interpolated terms denoted by y. 
For further details, the reader should consult the original work, in the 
Me'moires de PAoade'mie de Saint Pe'tersbourg, 1859. 

*The observations were taken in the years 1838 — 40, and tlie results were published in a 
long series of Uwaisons, both of the text and the charts, at different times for many years 
after the return of the expedition. 

fin case observations are equidistant and very near together, each area is regarded as 
equal to the sum of all the observations falling within its limits, multiplied by the constant 
interval between two consecutive observations. 



